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Polynomial parametrization of Pythagorean quadruples, quintuples and sextuples. 
Sophie Frisch (TU Graz) and Leonid Vaserstein (PSU) 



Abstract. For n = 4 or 6, the Pythagorean n-tuples admit a parametrization by a single n-tuple of 
polynomials with integer coefficients (which is impossible for n — 3). For n = 5, there is an integer- valued 
polynomial Pythagorean 5-tuple which parametrizes Pythagorean 5-tuples (similar to the case n = 3). 
Pythagorean quadruples are closely related to (integer) Descartes quadruples, which we also parametrize by 
a Descartes quadruple of polynomials with integer coefficients. 



Introduction 

An Pythagorean triple is a triple of integers (aq, x 2 , x 3 ) satisfying More generally, for any 

integer n > 3, and any commutative ring A, a Pythagorean n-tuple over A is an n-tuple (aq, . . . , x n ) 6 A n 
such that 

x 1 + • • • + x n _ 1 = x n . (1) 

Whenever A is not specified, we will understand A = Z. Likewise, "polynomial Pythagorean n-tuple" means 
a Pythagorean n-tuple over a ring of polynomials in finitely many indeterminates with coefficients in Z. 

Instead of studying (1) directly, it is often convenient to substitute u = x n + x n -i, v = x n ~ x n-i, and 
to consider the equation 

EiH r-a£_ 2 =w. (2) 

Over any ring A in which 2 is not a zero-divisor, this substitution and the reverse substitution 

x n -i = (u-v)/2, x n = (u + v)/2, (3) 

establish a bijection between solutions (aq, . . . , x n ) £ A n of (1) and solutions (aq, . . . , x„_2, u, v) G A n with 
u — v £ 2A of (2). 

We recall the existing polynomial parametrizations of integer Pythagorean triples. It is well known that 
up to permutation of aq and x%, every Pythagorean triple has the form 

(xi,x 2 ,x s ) =yo{2yiy 2 ,yl-yl,yl + yl) (4) 

with Hi G Z. In other words, the set of integer Pythagorean triples is the union of /i(Z 3 ) and where 

h(yo,yi,y2) = yo(2yiy 2 ,yl -yl,yl + yl) 

and 

f2{yo,yi,y2) = yo(yf - y 2 2,2 yi y 2 ,yl + yl) 

are two Pythagorean triples over the polynomial ring Z[yo, J/i, 2/2]- We say that all Pythagorean triples are 
covered by two polynomial Pythagorean triples (in 3 parameters each). 
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It is easy to see that the intersection of /i(Z 3 ) and /2(Z 3 ) contains only the zero triple. We know that it 
is not possible to cover all Pythagorean triples by any one Pythagorean triple over Z[j/i, . . . , y m ] for any m 
[FV]. It is, however, possible, to cover all integer Pythagorean triples by a single Pythagorean triple over the 
ring of integer- valued polynomials in 4 indeterminates [FV]. An integer- valued polynomial is a polynomial 
with rational coefficients which takes integer values whenever the variables take integer values. 

The primitive Pythagorean triples (xi, X2, X3) with positive X3, are, up to switching x\ and X2, given by 
(4) with primitive (t/i, 2/2) € ^ 2 such that 2/1 + yi is odd. The set of such pairs (yi, 2/2) admits a polynomial 
paramctrization [V]. Thus, all primitive Pythagorean triples can be covered by 4 polynomial triples (in 95 
parameters each, see [V], Example 14). 

All positive Pythagorean triples are, up to switching of x\ and x 2 , given by (4) with integers 2/1 > 2/2 > 0, 
2/o > 0. The set of such pairs admits a polynomial parametrization using the fact that every positive integer 
can be written as a sum of 4 squares plus 1. Thus, the positive Pythagorean triples can be covered by 2 
polynomial Pythagorean triples in 12 parameters. 

It is unknown whether the set of positive primitive Pythagorean triples can be parametrized by a finite 
set of polynomial Pythagorean triples. 

1. Quadruples 

After the short discussion of Pythagorean triples in the introduction, we now address the case n = 4, in 
other words, Pythagorean quadruples. 

Pythagorean quadruples were described by Carmichael [C], Chpt. II, §10, as follows: up to permutation 
of x\, X2, X3, every Pythagorean quadruple has the form 



with integer values for the parameters yo, . . . , 2/4. Thus, all Pythagorean quadruples are covered by 6 poly- 
nomial Pythagorean quadruples (in 5 parameters). Considering the position of the odd entry, it is easy to 
see that at least 3 permutations of xi,X2,X3 are needed. If one examines Carmichael's proof, one sees that 
three polynomial quadruples suffice, namely (xi, x 2 , x 3 , X4), (xi, X3, x 2 , X4) and (x 3 , x 2 , Xi, X4). 

We now show that a single polynomial Pythagorean quadruple covers all Pythagorean quadruples. Our 
proof does not make use of Carmichael's result (but rather provides a shorter proof of Carmichael's result 
as a byproduct). Nor do we use unique factorization in the ring of Gaussian integers Z[i] (which could be 
used to give an alternative proof). 

Definition. An n-tuplc w = (w\, . . . , w n ) S A n is called unimodular if w\A + ■ ■ ■ + w n A = A. In the case 
when A = Z this means that gcd(wi, . . . , w n ) — 1, i.e., w is primitive. 

Proposition 1. The integer solutions of 



(xi,x 2 ,x 3 ,x 4 ) = f(y , 2/1, 2/2, 2/3, 2/4) = 

2/0(22/12/3 + 22/22/4, 22/12/4 - 22/22/3, 2/1 + 2/2 - 2/3 - 2/4, 2/1 + vl + vi + yi) 



(5) 



2 1 2 
x l + x 2 = uv 



(6) 



are parametrized by the polynomial quadruple 



(x 1 ,x 2 ,u,v) = 2/0(2/12/3 + 2/22/4, 2/12/4 - 2/22/3, 2/1+2/2, 2/3+2/4) 



(7) 



as the parameters vary through the integers. Also, y can be restriced to odd integers. 
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Proof. We represent the integer solutions of (6) as Hcrmitian matrices 

_ / U X\ + 1X2 \ _ ( u X 

\X\ — 1X2 V J \X V 

of determinant over the Gaussian integers Z[z]. 

The group GL(2,Z[i]) acts on the Hermitian matrices as follows 

w -> g*wg (8) 

where g £ GL(2,Z[i]) and * means transposition composed with entry-wise action of complex conjugation. 
In particular, for an elementary matix g = E\2 (A) with A = Ai + A 2 i 



/l OUii iUl A\ / u (xi + Aim) + (x 2 + X 2 u)i \ 

\X w/V° 1 J ~ + X lU ) - (x 2 + X 2 u)i v + (Xj + \%)u + 2(\ 1 xi + X 2 x 2 ) J 

Setting either Ai = or A 2 = 0, we see that we can add an arbitrary integer multiple of u to x 2 , leaving 
u and x\ unchanged, and we can add an arbitrary multiple of u to n, leaving u and X2 unchanged. 

Given any solution with ti ^ 0, we can, by applying elementary matrices g in GL(2, Z[i]), make xi |, \x2\ < 
\u\/2, and hence \v\ < \u\ /2. Using the nontrivial permutation matrix in GL{2, Z[z]), we can switch u and v. 

Therefore, by an argument of descent, the orbit under GL(2,Z[i]) of any unimodular solution w = 
(x\,X2,u, v) to (6) contains a solution with v = (and hence x\ = X2 = 0) and u = 1 or —1. So we get 

with c = ±1, where (a, b) is the first row of the matrix g £ GL(2, Z[i]). 
So every integer solution of (6) has the form 

/ u xi + x 2 i\ ( a\ . . 
\Xi - x 2 i v J \b J 

with a, b £ Z[i], c G Z. Conversely, every expression of this form is a solution to (6) - it is not necessary to 
restrict (a, b) to be primitive or c to be ±1 or (sum of 2 squares)-free. 

Writing c = y , a = y\ + y 2 i and b = y^ + y±i, with indeterminates yk, we obtain a polynomial solution 

(U,V,X 1 ,X 2 ) = 2/o(2/l +Vl, Vl + Vl, 2/12/3 + ^2^4, 2/12/4-2/22/3) 

to (6) which covers all integer solutions. If we replace (a, b) above by (1 + i)(a, b), the solution is multiplied 
by 2. We can, therefore, restrict 2/0 to odd integers. □ 

Theorem 1. Let 

f(yo, 2/1,2/2,2/3, 2/4) = 2/0(22/12/3 + 22/22/4, 22/12/4 - 22/22/3, 2/? + yl - 2/1 - 2/4, S/i + V% + vi + vl) 
The polynomial Pythagorean quadruple 

.9 = /(2/o/2, 2/i, 2/2, 2/3, 2/1 +2/2 +2/3 + 2z) £ Z [2/0, 2/1, 2/2, 2/3, 2] 

in 5 parameters covers all Pythagorean quadruples, i.e., the range of the function g: Z 5 — > Z 4 consists of all 
Pythagorean quadruples. 
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Proof. Applying (3) to Proposition 1, we see that every Pythagorean quadruple has the form 

{x 1 ,x 2 ,x 3 ,x 4 ) = y (y 1 y 3 + y 2 y4, ym -2/22/3, (2/1 + vl - 2/3 - 2/I)/2, {y{ + y\ + y\ + 2/I)/ 2 ) 

where and 2/0(2/1+2/2+2/3+2/4) is even. Since y can be chosen odd, we may assume that 2/1 + 2/2 + 2/3 + 2/4 

is even. Writing y 4 = y\ + y 2 + 2/3 + 2z, we prove Theorem 1. □ 

To get Carmichael's result, note that x 4 is odd for any primitive Pythagorean quadruple (xi, x 2 , X3, X4) 
and that exactly one of xi,x 2 ,x 3 is also odd. So we can make x 3 + X4 even by switching, if necessary, X3 
with x\ or x 2 . Then gcd(xi, x 2l u, v) — 2 for the corresponding solution (xi, x 2l u, v) of (6). Going back from 
(7) with 2/0 = 2 to the Pythagorean quadruple, we obtain Carmichael's formulas. 

Notice that these formulas with yo = I and primitive (yi, y 2l 2/3, y 4 ) do not necessary give primitive 
solutions. Our proof shows that the necessary and sufficient condition for primitivity is the primitivity of 
(a, b) = (yi + y 2 i, y 3 + y 4 i). The set of primitive pairs of Gaussian integers admits a polynomial parametriza- 
tion by methods of [V], but this is beyond the scope of the present paper. 



2. Sextuples 



We discuss Pythagorean sextuples before quintuples because we will use sextuples in the proof of the 
parametrization of Pythagorean quintuples by a single quintuple of integer- valued polynomials in the next 
section. 

Dickson [Di], Section 106, attempted to describe all Pythagorean sextuples, i.e., all integer solutions to 

x\ + . . . + x\ = x\ (10) 

He observed that every integer solution of (10) gives rise to an integer solution of 

x\ + . . . + x\ = uv, (11) 

He solves the equation (11), in a lengthy proof of some 6 pages, but then fails to address the question of the 
reverse substition: how to return to Pythagorean sextuples from integer solutions of (11). 

We will also start by parametrizing the integer solutions of (11), giving a short proof using quaternions. 

Definition. The algebra of Lipschitz quaternions is the Z-algebra L generated by two symbols i and j 
subject to the defining relations i 2 = —1, j 2 = —1, and ji = —ij. We set k = ij. 

We recall a few facts about the algebra of Lipschitz quaternions. L is a free Z-module with basis 1, i, j, k 
and a free Z[i]-module with basis l,j. L can be represented as an algebra of 4 x 4 integer matrices or as an 
algebra of 2 x 2 matrices over Z[i] by identifying w = a + bi + cj + dk with 



M 4 (w) = 



( 


a 


b 


c 


d 


-b 


a 


-d 


c 




—c 


d 


a 


-b 


\ 


-d 


—c 


b 


a 



or M 2 (w) 



'! ■ hi c + di 
-c + di a — bi 



respectively. 

An involution on L is given by the Z-algebra anti-isomorphism 

a + bi + cj + dk (a + bi + cj + dk)* = a — bi — cj — dk. 

In the 4x4 integer matrix representation this corresponds to transposition; and in the 2x2 Gaussian integer 
matrix representation, to transposition followed by complex conjugation. 
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Definition. The norm of w = a + bi + cj + dk G L is defined as 

(a 2 + b 2 + c 2 + d 2 ) 2 = det(M 4 (w)) = (det M 2 {w)) 2 

and the reduced norm as 

a 2 + b 2 + c 2 + d 2 = w*w = det M 2 (w). 

For a 2 x 2 matrix M over L we define the norm of M as the determinant of the 8x8 integer matrix 
obtained by replacing each matrix entry w by M^w), and the reduced norm as the determinant of the 4x4 
matrix over Z[i] obtained by replacing each matrix entry w by M 2 (w). 

Remark. If w is a Hcrmitian 2x2 matrix over L, its entries commute and we can calculate the determinant 
in a naive way, as 

detf . u . , Xl + X2l + x ^ + x * k )=uv- x 2 -4-4~*i 

\xi — x 2 i — x 3 j — x 4 k v J ° 

The reduced norm of w as defined above is the square of this determinant. 
Proposition 2. A parameterization of all integer solutions (xi,x 2 ,x 3 ,X4,u, v) of 

x\ + . . . + x\ = uv 

in 9 parameters is given by 

xi = yo(yiV5 + 2/22/6 + 2/32/7 + 2/42/s) 
X2 = 2/0 (-2/12/6 + 2/22/5 + 2/32/8 - 2/42/7) 

xs = 2/0 (-2/12/7 - 2/22/8 + 2/32/5 + 2/42/e) 
Xi = 2/o(-2/i2/8 + 2/22/7 - 2/32/6 + 2/42/5) 

u = 2/o(2/? + 2/1 + 2/3 +2/4) 

v = Voivl +2/6 +V7 +2/1), 

as the parameters yo, ■ ■ ■ , J/8 range through the integers. Here yo may be restriced to ±1. 

Proof. We identify integer solutions w = (xi,x 2 ,x 3 ,X4,u,v) of x\ + . . . + x\ = uv with 2x2 Hcrmitian 
matrices over the algebra of Lipschitz quaternions L 

_ ( u xi + x 2 i + x 3 j + x 4 k \ 

\ xi — x 2 i — x 3 j — Xik v J 

of reduced norm 0. 

The group GL(2, L) acts on the set of 2 x 2 Hcrmitian matrices over L of reduced norm by (g, w) h-> g*wg, 
for g G GL(2,L), where g* results from g by application of the involution * to each entry, followed by 
transposition. 

Given any unimodular solution of (11) with u ^ 0, using an elementary matrix in GL(2, L), we can make 
\xi\, \x 2 \, \x 3 \, \x 4 \ < \u\/2, and hence \v\ < \u\. The inequality is strict unless \x m \ — \v\/2 for m = 1,2,3,4, 
in which case \v\ = 2 by unimodularity. In this last case, using an elementary matrix, we can arrange x m = 1 
for m = 1,2,3,4. 

Using the nontrivial permutation matrix in GL(2 1 L) 1 we can switch u and v. Therefore, by induction 
on |u|, the orbit of any unimodular solution w = (xi,x 2 ,x 3 ,X4,u,v) to (11) contains a solution with either 
\u\ = 1 and x m = for m = 1, 2, 3,4 or |u| = 2 and x m = 1 for m = 1, 2, 3, 4. 
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So we get that either 

with c = ±1 where (a, b) is the first row of the matrix g £ GL(2, L) or 

x * ( 2 l + i+j + k\ , /V \ , M 

™ = ±9 [ 1 _ i _ j _ k 2 ) 9 = ±^)(a,b) 

where (a, b) = ±(1 — i, 1 + j)g with g £ GL(2, L), because 

( 2 1 + i + j + k\ . • 1 , -wi ■ -i i -\ 

U-i-j-fc 2 J=(l-*,1+J) 

So every integer solution of (11) has the form 

/ u xi+ x 2 i + x 3 j + x 4 k \ ( a*\ , , . 

y xi — x 2 i — xsj — x 4 k v J \" J 

with a,b £ L and c £ Z. Here we need not restrict (a, b) to be primitive. 

Writing c = yo,a = yi+yii+Vsj + yik and b = 2/5 + 2/6* + 2/7 j + w e obtain the desired parameterization 
of all solutions w = (x\,X2,X3,X4,u,v) of (11). 

If we replace (a, b) above by d(a, b) with d £ L, the solution w is multiplied by d*d, which is equivalent to 
replacing y by yod*d. Since every nonnegative integer is of the form d*d (sum of 4 squares) we can restrict 
2/o to be ±1. □ 

Returning to Pythagorean n-tuples, the following polynomial Pythagorean sextuple is known: 

(xi,...,x 6 ) = h(y , . . . ,y 8 ) £ Z[y , . . . , y$} 6 (12) 

with 

xi = 22/o(2/l2/5 + 2/22/6 + 2/32/7 + 2/42/s) 
X2 = 2y Q (-y 1 y 6 + 2/22/5 + 2/32/8 - 2/42/7) 
X3 = 22/o(-2/i2/7 - 2/22/8 + 2/32/5 + 2/42/e) 
x 4 = 2y (-y 1 y 8 + 2/22/7 - 2/32/6 + 2/42/4) 

^5 = 2/0(2/? + 2/2 + 2/| + 2/4 - _ vl ~ V7 ~ vl) 
x e = yo(yi + 2/1 + vl + 2/4 + vl + 2/6 + Vi + vl) 

There are, however, integer Pythagorean sextuples such as (1, 1, 1, 1, 0, 2) that do not arise from the above 
polynomial sextuple with integer parameters 7/i- We now give a parametrization of all integer Pythagorean 
sextuples by a single polynomial Pythagorean sextuple in 9 parameters, or, by restricing the parameter yo 
to ±1, a parametrization by two integer Pythagorean sextuples in 8 parameters each. 

Theorem 2. Let h = h(y 0} ■ . ■ , yg) £ Z[y , . . . ,2/s] 6 be the polynomial Pythagorean sextuple (12) above. 
Then the polynomial Pythagorean sextuple in Z[yo, ■ ■ ■ , 2/7, z} 6 

5(2/o,2/i7 2/2,2/3,2/4,2/5,2/6,2/7,2) = /i(2/o/2, 2/1, 2/2, 2/3, 2/4, 2/5, 2/6, 2/7, 2/i + 2/2 + 2/3 + 2/4 + 2/5 + 2/6 + 2/7 + 2z) 

in 9 parameters covers all Pythagorean sextuples, i.e., the range of the function g: 1? — > Z 6 is precisely the 
set of all Pythagorean sextuples. Also, the parameter y can be restricted toil. 
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Proof. We obtain all Pythagorean sextuples from all solutions of (11) with u — v even by by (3). Since 
we we can take j/o odd (even ±1) in Proposition 2, we may assume that y\ + • • • + yg is even. Writing 
2/8 = 2/1 + • • • + 2/7 + 2z, we obtain a Pythagorean sextuple over 1\yo, • • ■ , 2/7, z] in 9 parameters which 
parametrizes all Pythagorean sextuples: 

xi = 2/o(?/i2/5 + 2/22/6 + 2/32/7 + 2/4(2/1 H h 2/7 + 2z)) (13) 

£2 = 2/0 (-2/12/6 + 2/22/5 + 2/3(2/1 H 1- 2/7 + 2z) - 2/42/7) 

£3 = 2/o(-2/i2/7 - 2/2(2/1 H h 2/7 + 2z) + J/3J/5 + y 4 y 6 ) 

2:4 = 2/o(-2/i(2/i H 1- 2/7 + 2z) + 2/22/7 - 2/32/6 + 2/42/5) 

x 5 = 2/0(2/? + 2/ 2 2 + 2/ 3 2 + 2/4 - 2/5 - yl - y 2 7 - (2/1 + • • • + 2/7 + 2z) 2 )/2 
^6 = 2/o(2/i + 2/ 2 2 + 2/ 3 2 + 2/4 + 2/5 + vl + y 2 7 + (2/1 + • • • + 2/7 + 2z) 2 )/2 

□ 

Remark. The Z- algebra L is a subring of the ring H of rational quaternions (or Hamilton quaternions). 
Adjoining to L the element + + k)/2, we obtain a larger subring L' of H, called the ring of Hurwitz 
quaternions. This ring H has certain unique factorization properties, which, however, we did not use. 
They could be used to give alternative proofs for the results in this section. 

3. Quintuples 

We now consider the case n = 5 of Pythagorean quintuples. We obtain Theorem 3 from Proposition 2 
via Proposition 3. 

Proposition 3. A parametrization of all integer quintuples (x\,X2,Xs, u, v) satisfying 

2 , 2 , 2 

x 1 + x 2 + x 3 = uv 

by a quintuple of polynomials with integer coefficients in the 12 parameters yo, zq, z\, Z2, Z3, Z4, Z12, 213, 
214, ^23, z 2 4, Z34 is given by 

x\ = 2/0(2/12/5 + 2/22/6 + 2/32/7 + 2/42/8) 
X2 = 2/0 (-2/12/6 + 2/22/5 + 2/32/8 - 2/42/7) 
X3 = 2/0 (-2/12/7 - 2/22/8 + 2/32/5 + 2/42/6) 

u = 2/0(2/1 +2/2 +2/3 +2/4) 

« = 2/o(2/ 2 +2/ 2 + 2/t ! +2/8)> 

with 

yi=z z 1 , y2 = z z 2 , 2/3 = z z 3 , y 4 = z z 47 

2/5 = -ZuZ\ - Z 2 aZ2 - Z34Z3 
2/6 = 213,21 + Z2322 - 234Z4 
2/7 = -21221 + 22323 + 22424 
2/8 = -21222 - 2i 3 2 3 - 21424. 
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Proof. To solve x\ + x\ + x\ = uv, we set x 4 = in the general solution to x\ + x\ + x\ + x\ = uv obtained 
in Proposition 2. 

-2/12/8 + 2/22/7 - 2/32/6 + 2/42/5 = (14) 

(The case ?/o = only contributes the zero solution which we will not miss.) 

The integer solutions of (14) can be parametrized by 11 parameters as follows. 

First we write (2/1,2/2,2/372/4) — 20(21, z 2, 2 3 , 24) with £ Z and (21,22,23,24) unimodular such that 
(21, Z2, 23, 24). (-2/8, 2/7, -J/6,2/5) = (the case 2/1=2/2 = 2/3 = 2/4 = only contributes solutions (0,0, 0,0, v) 
which we will not miss). 

By [VS], Remark after Lemma 9.6, we can write 

(-2/8,2/7, -2/6,2/5) = 

212(2:2, -21, 0, 0)+2l 3 (2 3 , 0, -2l, 0)+2l 4 (2 4 , 0, 0, -2l)+2 23 (0, 2 3 , -2 2 , 0) + 2 24 (0, 2 4 , 0, -2 2 )+2 34 (0, 0, 2 4 , -2 3 ). 

This gives a parametrization of the integer solutions of (14) in the 11 parameters 2 , 21, 2 2 , 23, 2 4 , 212, 

213, 2l 4 , 2 23 , 2 24 , 2 34 . 

Therefore all integer solutions of x\ + x\ + x\ = uv are parametrized by a polynomial solution with 12 
parameters including 2/0- □ 

Another parametrization of x\ + x\ + x\ = uv with 20 parameters can be obtained using [V] , Proposition 
3.4 with k — 8. We now parametrize integer Pythagorean quintuples by a single Pythagorean quintuple 
over the ring of integer- valued polynomials in 14 variables. This can be used to construct a parametrization 
by a finite number of integer-coefficient polynomial Pythagorean quintuples [F]. Whether it is possible to 
parametrize integer Pythagorean quintuples by a single quintuple of integer-coefficient polynomials or not, 
we do not know. 

Theorem 3. A parametrization of all Pythagorean quintuples by a quintuple of integer-valued polynomials 
in the 14 variables w , w 12 , w 13 ,w 14 , w 23 , w 2 4, w 34 , h,t 2 , t 3 , d 1} d 2} d 3 , w 4 is given by (/1, / 2 , / 3 , f 5 , f 6 ), where 

h = 22/o(2/i2/5 + 2/22/6 + 2/32/7 + 2/42/8) 
h = 22/o(-2/i2/6 + 2/22/5 + 2/32/8 - 2/42/7) 
h = 22/o(-2/i2/7 - 2/22/8 + 2/32/5 + 2/42/e) 

f I 1 1 2 , 2 , 2 2 2 2 2\ In 

h = 2/o(2/i + 2/ 2 + 2/3 + 2/4 - 2/5 - 2/ 6 - 2/7 - 2/ 8 )/2 
/a = 2/o(2/? + 2/ 2 + 2/ 3 2 + 2/1 + 2y| + Vl + 2/? + 2/D/2 

and 

2/l= 2 2i, 2/2 = 2 2 2 , 2/3 = 2 2 3 , 2/4 = 2 2 4 , 
2/5 = -21421 - 2242 2 - 2 34 2 3 
2/6 = 2l 3 2l + 2 23 2 2 - 2 34 2 4 
2/7 = -21221 + 2 23 2 3 + 2 24 2 4 
2/8 = -2l 2 2 2 - 2 i3 2 3 - 2 i4 2 4 . 

and further 
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z =w + tiw + t 2 w - 2tit 2 w + t 3 w - 2tit 3 w - t 2 t 3 w + 2tit 2 t 3 WQ + txw 12 - tit 2 w 12 - 

- tihw\ 2 + t 2 t 3 wi 2 + t 2 wi 3 - ht 2 wi 3 + t 3 wu - tit 3 wi4 + tiw 23 + t 2 w 23 - 2tit 2 w 23 - 

- ht 3 w 2 3 - t 2 t 3 w 23 + 2tit 2 t 3 w 23 + tiw 2i - ht 2 w 24 + t 3 w 2i - 2tit 3 w 24 , - t 2 t 3 w 24 + 
+ 2t 1 t 2 t 3 w 24 + t 2 w 34 - ht 2 w 3i + t 3 w 34 - ht 3 w 3i - 2t 2 t 3 w 34 + 2tit 2 t 3 w 3i 

zi =2di + t\t 2 + t 3 - 2ht 2 t 3 + Wi 
z 2 =2d 2 + ti - t x t 2 + t 3 - t x t 3 - t 2 t 3 + 2t 1 t 2 t 3 + w 4 
z 3 =2d 3 +t 2 + t 3 - ht 3 - 2t 2 t 3 + 2ht 2 t 3 + w 4 
z 4 =w 4 

Z12 =wi 2 + ht 2 wi 2 - ht 2 t 3 w 12 + ht 2 wi 4 - ht 2 t 3 wi 4 + ht 2 w 23 - ht 2 t 3 w 23 + ht 2 w 34 - ht 2 t 3 w 34 

Z13 =W13 + t\t 3 Wl 3 - tit 2 t 3 W 13 + tlt 3 Wli - t\t 2 t 3 Wl 4 + t\t 3 W 23 - tit 2 t 3 W 23 + t\t 3 W 24 - t\t 2 t 3 W 24 
Z 14 =W!4 
Z23 =W 23 

Z24 =tit 2 t3Wi 2 + tit 2 t3W\3 + W 24 + t\t 2 t3W 24 + tit 2 t3W3 4 

Z34 =W 34 

Proof. To go from the solutions of x\ + x\ + x\ = uv parametrized in Proposition 3 to the solutions of 
x\ + x\ + x\ + x\ — x\ we use (3), allowing only those u, v with u ± v even. 

In our case, we need to parametrize those zo, . . . , z 34 that make yi+y 2 + - ■ -+y% even, i.e., those z , . . . , z 34 

such that E — z^z\ + z$z 2 + zqz 3 + z$z 4 H z\ 4 zi — z 24 z 2 — Z34Z3 — Z13Z1 — z 2 3Z 2 + Z34Z4 — zi 2 z\ + Z23Z3 + 

Z24Z4 - z 12 z 2 — Z13Z3 — Z14Z4 is even. 

This is acheived by the following substitution, which, after simplification, gives the parametrization in 
the statement of the theorem. 

(z , Zl, Z 2 , Z 3 , Z 4 , Z12, Z13, Z U , Z23, Z 24 , Z34) = 

(w ,w 4 + 2d ll w 4 + 2d 2 ,w 4 + 2d3,w 4 ,w 12 ,wi3,wi 4 ,w 2 3 1 w 24 ,W3 4 )(l - - t 2 )(l - t 3 )+ 
(w\i+w 2 i+w 3i +2wo,w i +2d\ + l,w i +2d 2 +l,w i +2d 3 + l, Wi,w\ 2 ,wi 3 ,w\i, w 23 ,w 2i , w 3 i){l-ti){l-t 2 )t 3 + 
(■W13 + w 23 + w 34 + 2w ,w 4 + 2dx,w 4 + 2d 2 , w 4 + 2d 3 + 1, w 4l w 12 , wi 3 ,w 14l w 23l w 24 , W34) (1 - ii)t 2 (l - t 3 )+ 
(wi 2 + w 23 + w 24 + 2w ,w 4 + 2di,w 4 + 2d 2 + l,w 4 + 2d 3 , w 4 , wi 2 , W13, w\ 4 , w 23l w 24 , w 34 )ii(l - t 2 )(l - t 3 ) + 

(wi2 + wi 3 + W14 + 2w , Wi + 2di + 1, Wi + 2d 2 ,w 4 + 2d 3l w 4l Wi 2 ,wi3,wi 4l w 2 3,w 24 , W34) (1 - h)t 2 t 3 + 
(w , w 4 + 2d\ + 1, w 4 + 2d 2 + l,w 4 + 2d 3 ,w 4l wi 2 ,w 23 + w 24 + w 14 + 2w 13 , w 14 , w 23l w 24 , ^34)^(1 - t 2 )t 3 + 
(w , w 4 + 2di + 1, u>i + 2d 2 ,w 4 + 2d 3 + 1, w 4 , w 23 + w 14 + w 34 + 2w i2 , W13, Wu, w 23l w 24 ,W3 4 )tit 2 (l - t 3 ) + 
(wq,w 4 + 2di,w 4 + 2d 2 + 1, Wi + 2d 3 + 1, w 4l wi 2 , wi 3l wu, w 23 , w 12 + w i3 + w 34 + 2w 24 ,W3 4 )tit 2 t 3 

□ 

4. Descartes quadruples 

In 1643 Descartes [De] described a relationship between the radii of four mutually tangent circles (called 
a Descartes configuration), namely, 

2{b\ + bj + b 2 3 + bl) = (61 + b 2 + b 3 + b 4 ) 2 (15) 
9 



where bi are the reciprocals of the radii. Others, including Steiner, Beecroft, and Soddy [S], rediscovered the 
result. We call an integer solution of (15) a Descartes quadruple. 

Given one Descartes configuration, there is a geometric way to produce plenty of them creating an 
Apollonian packing. If the four intitial curvatures bi are integers, all curvatures in the packing are integers. 
There are several publications about integer Apollonian packings [GLMWY1], [GLMWY2], [GLMWY3], 
[EL] , [N] , [LMW] . A bijection between integer Pythagorean quadruples and integer Descartes quadruples can 
be found in [GLMWY1], Lemma 2.1. 

In this section we parametrize all integer solutions of (15) by a single polynomial solution in 5 parameters, 
using a bijection between integer Descartes quadruples and integer solutions of (6). 

Given an integer solution (xi,x 2 ,u,v) of (6), 

bi = u + v — 2x\ + x 2 , b 2 = u + x 2 , b 3 = v + x 2 , & 4 = -x 2 (16) 

is an integer solution of (15). Conversely, we can invert this linear transformation: given an integer solution 
(bi, b 2 , &3, bi) of (15), b\ + b 2 + b 3 + 64 is even and 

x\ = (-61 + b 2 + H + bi)/2, x 2 = -bi, u = b 2 + b i , v = b 3 + b 4 

is an integer solution of (6). 

Theorem 4. A parametrization of all integer solutions of 

2(b\ + b 2 2 + b\ + b\) = {h + b 2 + b 3 + 6 4 ) 2 (15) 

in 5 parameters is given by 

b\ = 2/o(?/i + 2/1 + vl + Va - 2 2/i2/3 - 22/22/4 + 2/12/4 - 2/22/3), 
b 2 = 2/0(2/1 + 2/1 + 2/12/4 - 2/22/3), 

63 = 2/0(2/3 + 2/4 + vwi - 2/22/3), 

bi = 2/0 (-2/12/4 - 2/22/3)- 

Proof. In the expression (16) of &i, b 2 , 63, &4 in terms of a solution x\, x 2 , u, v of (6) we have substituted the 
parametrization of all integer solutions of (6) from Proposition 1. □ 
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